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We propose versatile materials based on the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 
in a zigzag silicene nanoribbon (ZSNR) on half filling in the presence of an out-of-plane electric field. 

We show that the topological phase transition in the band dispersion of ZSNR can be probed by 
using the RKKY interaction. We find that, due to the zero-energy edge states of the ZSNR, the 
exchange coupling is significantly enhanced when the impurities are located on the zigzag edges, and 
also explore that the strength of the interaction in the topological insulator phase is much greater 
than that when the system is in the band insulator region. We present a model to investigate the 
phase of a system of two magnetic impurities located on the edge of the ZSNR and find that three 
different magnetic phases, spiral, ferromagnetic, and anti-ferromagnetic, are possible for different 
values of the electric field. This electrical tunability of the magnetic phases in silicene can be 
explored by using current experimental techniques and can be of interest in the field of spintronics. 

PACS numbers: 68.65.Pq, 73.23.-b, 81.05.Zx 


I. INTRODUCTION 

The idea of spintronics is to utilize the spin degree of 
freedom of an electron associated with its intrinsic an¬ 
gular momentum [1, 2]. This field is expected to play 
an important role in future technologies, such as infor¬ 
mation storage, processing and communication at faster 
speed and lower-energy consumption [3] . Among the dif¬ 
ferent candidates in the field of spintronics including the 
d and / blocks of the periodic table, carbon and silicon 
have their own attractions, such as low density, biocom¬ 
patibility, and plasticity. Although these atoms are not 
magnetic, there are some proposals to use their deriva¬ 
tives in the field of spintronics. One of the proposals is 
to use a doped system with magnetic adatoms. Magnetic 
impurities in dilute magnetic semiconductors can interact 
indirectly via the itinerant electrons of the host material, 
which is known as the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction [4]. This interaction is a vital inter¬ 
action and plays a significant role in the magnetic order¬ 
ing of many electronic systems, including spin glasses [5] 
and alloys [6]. 

The RKKY interaction has been extensively studied 
in carbon based materials [7-17] and it can be experi¬ 
mentally measured from single-atomic magnetometry of 
a pair of magnetic atoms [18, 19] and also magnetotrans¬ 
port measurements [20] based on the angle-resolved pho¬ 
toemission spectroscopy (ARPES) and scanning tunnel¬ 
ing microscopy (STM) experiments. While this interac¬ 
tion has a simple Heisenberg term for spin nondegener¬ 
ate materials such as graphene, it contains an additional 
Ising term in spin polarized graphene [12]. The Ising and 
Heisenberg terms of the RKKY coupling make the mag¬ 
netic moment of impurities to be collinear. However, the 
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Rashba or Dresselhaus terms lead to a Dzyaloshinskii- 
Moriya (DM) term, which causes a rotation between the 
spinors of the magnetic impurities. This phenomenon 
was shown to appear in three-dimensional topological in¬ 
sulators [21-23] as well as a two-dimensional electron gas 
with the Rashba term [24]. It has been also shown that 
the DM term can appear in a monolayer of molybdenum 
disulfide [25] and a graphene nanoribbon together with 
spin-orbit coupling [17], where the Rashba term is no 
longer present and the spin-valley coupling is the reason 
behind the twisted interaction between the impurities. In 
addition, the existence of the DM interaction has been 
reported for a bulk silicene [26]. 

Unlike carbon atoms in graphene, silicon atoms in 
silicene tend to adopt the tetrahedral sp^ hybridiza¬ 
tion over sp^, which results in a slightly buckled struc¬ 
ture of silicene, where sublattices are displaced from 
each other in the out-of-plane direction with a distance 
Takeda, Durgun, Cahangirov. This buckled structure 
also leads to a relatively large spin-orbit coupling (SOC) 
{^so = 3.9 meV) [30] which is about 1000 times larger 
than graphene [31, 32]. The buckled honeycomb lattice 
and its Dirac band structure, combined with its sensitive 
surface, offers the potential for a widely tunable two- 
dimensional monolayer, where external fields can be ex¬ 
ploited to influence the fundamental properties for future 
nanoelectronic devices [33, 34]. The buckled structure of 
silicene makes it possible to control its band dispersion 
by using an external out-of-plane electric field. In this 
scheme, the on-site energies of different sublattices are 
no longer equivalent and thus the spin-up and spin-down 
bands become separated [35, 36]. 

A transition from a topological insulator state to a 
band insulator in silicene is predicted at E^dj^ = 3.9 
meV by increasing the strength of the electric potential 
[35, 37-39]. Although there exists a finite gap in the bulk 
dispersion in both regions, the bulk-edge correspondence 
leads to a zero-gap structure in the topological phase. 
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while such zero-energy edge states disappear in the band 
insulator state. Previously, the effect of this topological 
phase transition had been explored in silicene via Friedel 
oscillations related to the existence of a single impurity 
[40]. In this work, we are aiming to find the effect of 
this topological phase transition in the nonlocal interac¬ 
tion between two magnetic impurities. We show that the 
topological phase can be probed using the RKKY interac¬ 
tion. To do so, we consider a finite-size ZSNRzigzag (see 
Fig.l) and consider impurities to be located on the edge 
and calculate the RKKY interaction as a function of the 
distance between isilicene nanoribbon (ZSNR)mpurities 
for two different phases of normal and topological band 
insulators. To show that our findings are affected by 
the topological phase of the ZSNR, we compare our re¬ 
sults with those when the impurities are located in the 
bulk. Moreover, to propose a practical application, we 
assume a system of two magnetic impurities on the edge 
of ZSNR and demonstrate that the ground state of such 
a system can displays spiral, ferromagnetic (FM), or an¬ 
tiferromagnetic (AFM) phases depending on the value of 
the out-of-plane electric field. Each of these phases has 
its own attraction in the field of spintronics [41]. Because 
of the electrical tunability of the magnetic phases of the 
system which is due to the electronic band tunability of 
silicene, it is of much interest to use magnetically doped 
silecene in nanotechnology [42]. 

The paper is organized as follows. In Sec. II, we first 
show the tight-binding Hamiltonian of a silicene nanorib¬ 
bon in the presence of an out-of-plane electric field and 
then introduce the RKKY interaction in terms of non¬ 
interacting Green’s function of the ZSNR. In Sec. Ill, 
we present our numeric calculations for the coupling 
strengths of the RKKY interaction in the ZSNRs in terms 
of distance between the impurities. Also, different pos¬ 
sible phases of magnetism for a system of two magnetic 
impurities on the edge of silicene are discovered. Sec¬ 
tion IV contains discussions and a brief summary of our 
main results. 


II. THE MODEL HAMILTONIAN AND 
FORMALISM 

The low-energy excitation of silicene is described by 
Dirac electrons with a large SOC interaction. The phase 
diagrams of silicene as functions of the electric potential 
and exchange field exhibit a rich physics [43], for instance, 
a quantum anomalous Hall insulator, valley polarized 
metal, quantum spin Hall, and band insulator phases ap¬ 
pear. Remarkably, they are characterized by the Chern 
number and spin-Chern number together with the edge 
modes of a nanoribbon. The tight-binding Hamiltonian 
for a single layer silicene exposed to an external out-of¬ 
plane electric field Ez can be written as [30, 35] 


+ ;YCi^kV., (1) 

is 

where <> and run over all the nearest and next- 
nearest neighbor hoping sites, respectively. c\^ {ci^) cre¬ 
ates (annihilates) an electron with spin a at site i, and 
t = 1.6eV is the transfer energy. Uij = dij • z/\dij \ = ±1 
{dij = di X dj) where d^,dj are the vectors from center 
of a hexagon to the i and j adjacent to the same sublat¬ 
tices, respectively. = ±1 for the A or B site, a are the 
Pauli spin matrices and I = d/2. By performing a Fourier 
transformation in the continuum model, the low-energy 
effective Hamiltonian in the vicinity of the K or K' point 
is 

nik) = hvpikxTx - kyTyni) - XsoTzCTzil + IEzTz, ( 2 ) 

where tj = ±1 for the K or K' Dirac point and cr,, r* 
are Pauli matrices for spin and sublattice pseudospin, 
respectively. The Fermi velocity is up = with the 

lattice constant a = 3.86A. The energy spectrum is thus 
given by 

Er,{k) = ±^{nvFk)^ + {IE, - s.vXso?- ( 3 ) 

In the absence of an external electric field, there are two 
separated bands with a SOC gap of 2Ago- However, by 
applying Ez , the inversion symmetry breaking causes the 
spin splitting of \lEz ± t]SzXso\ [44], where Sz = El 
denots the spin of electrons. The gap energy Eg de¬ 
creases linearly with increasing Ez until a critical field 
Eq = T]S zXgo 11 - 

Since it is possible to engineer the electronic structure 
of two-dimensional (2D) nanorribons by using scanning 
tunneling microscope lithography with nanometer preci¬ 
sion of the patterning of nanoribbons [45], we consider 
a zigzag silicene nanoribbon [46, 47] with a finite width 
along the y— direction, where its properties can be de¬ 
scribed by Eq. (1). A schematic representation of the 
ZSNR is shown in Eig. 1. The dashed box is a unit cell 
and index m represents the number of unit cells or equiv¬ 
alently the X coordinate, and n shows the y coordinate 
of the lattice points, and therefore each site can be la¬ 
beled by (m, n), a notation that we use here to show the 
position of the impurity. 

By locating magnetic impurities on the ZSNR, an ex¬ 
change interaction between the spin of the magnetic im¬ 
purities and the spin of the conduction electrons of the 
silicene atoms occurs, which is given by 

^int = Jc Y (4) 

i=i,2 

where Jc is the coupling constant between the conduc¬ 
tion electrons and impurities, Ij is the magnetic mo¬ 
ment of the j-th impurity at location Rj, and S(R) = 
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and Ising terms, respectively, and all the other parts of 
the spin susceptibility are zero. The explicit form of Xij 
as a function of eigenvectors and eigenvalues of the sys¬ 
tem is given in the Appendix. 

III. NUMERICAL RESULTS 

The general features of the exchange coupling, basi¬ 
cally the dependence of the RKKY interaction on the 
distance R, have been numerically studied previously [26] 
for a bulk silicene system. Here, we mainly focus on the 
half filling (ef = 0 ) for various configurations for which 
the positions of two magnetic impurities are either in¬ 
side the ZSNR or on the edges. Our main purpose here 
FIG. 1. (Color online) Lattice structure of the zigzag silicene explore the topological phase and edge states de- 

nanoribbon. The edge and inside configurations marked ^ence of the RKKY interaction. In this regard, we 

With biack and green spins, respectively, ihe dashed box . , , i • p i i i i • 

., 11 ? , j j- . c consider two electric field strengths m which their corre- 

is a unit ceil and m^n represent the x and y coordinates ot Tin i i 

the lattice points. The length (L) and the width (A) of the ^ponding bulk gaps are equal. For instance, we choose 
nanoribbon are defined as the number of unit cells and the — 1 meV, and lEz — 6.8 meV, where the bulk band 

number of atoms in a unit cell, respectively. Here, we have gaps of these two fields are both 2.9 meV, and then cal- 

taken A = 12 . culate the RKKY interaction between the two magnetic 

impurities as a function of distance between them for dif¬ 
ferent impurity configurations. 

2 Si ~ where (Ji is the spin density operator Figures 2 (a) and 2 (b) show the spatial behavior of 

of the conduction electron. RKKY interaction for two impurities, both sitting on 

Using a second-order perturbation theory, one can re- same edge along the zigzag direction for A = 12 and 

place the Hint by the RKKY interaction given by [4, 24] ^ ^ 300 Fig. 2 (a)/A^ = 1 meV and Fig. 2 (b)/A^ = 6.8 

meV. Figures 2 (c) and 2 (d) are the same as Figs. 2 (a) 
and 2 (b), but for two impurities located inside the zigzag 
^ nanoribbon along the line n = 6 . As shown in this figure, 
Jrkky foi* all terms displays an oscillatory behavior with 
Here, a and /3 indicate the sublattices of which the respect to the distance between two impurities. As an 

impurities are located on and Xki {kj = x,y, and z) is important result, by comparing Figs. 2 (a) and 2 (c) 

the spin susceptibility tensor of the system in real space, and Figs. 2 (b) and 2 (d), the RKKY coupling in the 

obtained through the retarded Green’s functions cases shown in Figs. 2 (a) and 2 (b) is significantly larger 

than when both impurities are in the bulk due to the 
kA r 2 Q , f . I Q , f xn existence of nearly-zero-energy states at the edge of the 

^ tt J_^ cxp\ ^ w pa\ ^ ^ ^^’ZSNR [38]. This fact IS true even m the cases where the 

j^ 0 ^ystem is in the band insulator regime [Figs. 2 (b) and 2 
(d)] since the local density of states drops dramatically 
where sf is the Fermi energy, and the trace is over the by going from the edge to the bulk [49], and therefore 

spin degree of freedom, the In order to calculate the spin the RKKY interaction is much larger in the case shown 

susceptibility of the finite-size silicene nanoribbon, we use in Fig. 2 (b) in comparison with results demonstrated in 
the spectral representation of the noninteracting Green’s Fig. 2 (d). In addition, in the topological insulator region 
function [48], {lEz = 1 meV< Ago), the RKKY interaction is about 20 

times greater than in the band insulator region (lEz = 6.8 
meV> Ago) when impurities are located on the edge, as 
shown in Figs. 2 (a) and 2 (b), while no difference in the 
order of magnitude appears in the case where impurities 
where the sum runs over all eigenstates V’n.s and corre- are in the bulk, as shown in Figs. 2 (c) and 2 (d). Note 

spending eigenvalues with spin s and the band index that although Jrkky shows a few oscillations in R owing 

n obtained by diagonalizing the finite-size Hamiltonian, to the smallness of the gap of silicene, it decays fast with 

Eq. (1). The RKKY Hamiltonian, can be written as a short-ranged behavior. We also calculate the RKKY 

^ 7 TTi 7 \ z \ IT T(Q\ Strength in the case of one impurity on the edge (11,1) 

riRKKY = ^2) ' ZE JlllzEz^o) j r / 

and the other one on the inner edge on the line (m, 2 ) 
where Jh = J^Xxx^ Jdm = JcXxy and Jj = J^iXzz — with m = 12,13,... (gray arrow in Fig. 1 ). The obtained 

Xxx) represents the Heisenberg, Dzyalosinskii-Moriya, data are qualitatively similar to those results shown in 
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FIG. 2. (Color online) Various terms of the RKKY interaction 
scaled by as a function of the impurity distances for the 
ZSNR with N = 12, L = 300.(a),(b) Impurities located on 
the edge (black arrows) such that the hrst impurity is hxed at 
the sublattice (11,1) and the second one is located at (m,l) 
where m — 12,13,.... (c),( d) Impurities is located inside the 
ZSNR(green arrows) such that the hrst impurity located on 
the sublattice (11,6) and the second one is placed at (m, 6) 
with m=12,13,.... In these hgures, the electric held is hxed 
at lEz = 1 meV in (a) and (c) lEz = 6.8 meV in (b) and 
(d). Note that Jrkky shows an oscillatory behavior in R 
and decays fast with a short-ranged behavior. The RKKY 
coupling is signihcantly larger than the situation where both 
impurities are in the bulk due to the existence of nearly-zero- 
energy states at the edge of the ZSNR. 


Figs. 2 (c) and 2 (d) which explore that both impurities 
might be sitting on the edge to obtain the topological 
effect. 

We calculate the electric held dependence of the 
RKKY interaction for different distances between mag¬ 
netic impurities located on the edge to explore the effect 
of the topological phase transition on the RKKY inter¬ 
action. Figure 3 shows the electric potential dependence 
of the various terms of the RKKY coupling, when both 
impurities are located on the edge of ZSNR. To include 
the effect of different distances between impurities, we 
consider impurities with distances from la to 10a on the 
edge and take the root mean square value of the RKKY 
interaction. As seen in Fig. 3, in the topological insula¬ 
tor region {Ez < Ec), all the interaction coefhcients are 
much greater than those results obtained for the band 
insulator region {Ez > Ec). Note that this behavior is 
special to the edge of ZSNR and as one calculates the 
same for the bulk no change in the order of magnitude 
can be seen for different phases. 

We compare various terms of the RKKY interaction 
for two widths of the ZSNR in Fig. 4 to study the ef¬ 
fect of the width of ZSNR on the RKKY coupling. We 
show the RKKY interaction (scaled by J^) as a function 


xlO 



FIG. 3. (Golor online) Root mean square of various terms of 
the RKKY interaction scaled by as a function of the per¬ 
pendicular electric potential for the ZSNR with N — 22, L — 
300. Here, the hrst impurity is hxed at the sublattice (147,1) 
and the position of the second one is varied from sublattices 
(148,1) to (157,1). Ec in this hgure has been labeled by a 
black vertical arrow. 





FIG. 4. (Golor online) Various terms of the RKKY interaction 
scaled by as a function of the perpendicular electric poten¬ 
tial for L = 300, when two impurities are hxed at the zigzag 
edge at (148,1) and (153,1) lattice points (a) the Heisenberg 
and Ising terms and (b) the Dzyaloshinsky-Moriya term. In 
both cases R/a = 5. Noticeably, the RKKY coupling in¬ 
creases with ZSNR width due to the increase of the density 
of states of the edge mode. Black vertical arrows show the 
position of Ec. 
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of the perpendicular electric potential, for two impuri¬ 
ties located with a fixed distance on the zigzag edge at 
(148,1) and (153,1) lattice points. Figure 4 (a) shows 
the Heisenberg and the Ising terms and Fig. 4 (a) be¬ 
longs to the Dzyaloshinskii-Moriya term. In both cases 
R/a = b. By comparing the results of various systems 
with different widths, it is clearly found that the RKKY 
coupling increases with the ZSNR width, as reported pre¬ 
viously [38, 50]. This means that thicker nanoribbons 
can show better topological effects in ZSNRs. Although 
we assume a narrow nanoribbon in our numerical calcu¬ 
lations, we expect that sharper differences between the 
results of the bulk and edge states should be appear in 
the thicker cases too. In contrast to changing the width 
of ZSNR, we find that changing the length of the sample 
has little effect on the RKKY interaction. Note that since 
we assume the Fermi energy to be zero in these figures, 
all the diagrams are short ranged and fall off rapidly with 
an exponential like function. Most importantly, we ex¬ 
amine the RKKY interaction for an n-doped ZSNR and 
we find that the RKKY coupling behaves as R~^ in the 
finite values of the doping which shows that the edges of 
doped ZSNR are acting as a one-dimensional (ID) semi¬ 
conductor. 

Now, we assume a system of two magnetic impuri¬ 
ties located on the edge of ZSNR with constant distance. 
This would be a useful and practical application of the 
RKKY interaction in a nanoribbon of silicene. Our aim 
is to explore the final phase of such a system. Following 
Ref. [25], we assume magnetic impurities as classical 
spin vectors and rewrite the RKKY interaction Hamil¬ 
tonian in the spherical coordinate. We find that the 
equilibrium direction of the vectors of the magnetic mo¬ 
ments is 0i^2 = 0, TT or = $2 = 7r/2, where 0 is 
the polar angle and it is defined as an angle between 
the spin direction with the unit vector normal to the 
ZSNR surface (z). Also, we find that the impurity mo¬ 
ments shall be located inside the plane with a relative az¬ 
imuthal angle of 0 {Jdm/ Jh)^ if two quantities 

Di = —sign^Jn )’ ^2 = J^m~Ji) 
are both positive {Di and D 2 are quantities of the Hes¬ 
sian matrix introduced in Ref [25)]. Otherwise, the mag¬ 
netic moment of impurities will be perpendicular to the 
plane and will form a FM structure when Jh + J/ < 0 or 
an AFM structure for the case when Jh ^ 

Figure 5 shows the sign of Di and D 2 as a function 
of the out-of-plane electric potential lEz at half filling 
in the case where impurities are located on the edge of 
the nanoribbon with distance R = ba. A change of the 
equilibrium state of the system from the case where im¬ 
purities have a finite angle to the case where they have 
a parallel, FM-like structure, occurs at lEz ^ 0.91 meV. 
Afterwards, another equilibrium state occurs from paral¬ 
lel to antiparallel magnetic moments at lEz ^ 1.35 meV. 
It would be worthwhile mentioning that in the spiral re¬ 
gion, the direction of the magnetic moment is in the plane 
of the sample such that each impurity rotates with an an¬ 
gle (j) with respect to the next one. This angle of rotation 
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FIG. 5. (Color online) Sign (Di, 2 ) and phase diagram of a 
system of two magnetic impurities placed on the edge of ZSNR 
with N — 12, L = 300 as a function of the external potential. 
Here, impurities are located with five lattice vectors {R/a — 
5) from each other, at (148,1) and (153,1) lattice points. 


starts from zero and increases up to high values of about 
80A as Ez increases. 


IV. CONCLUSION 

In conclusion, we have studied the effect of the out-of- 
plane electric field and topological phase transition on the 
RKKY interaction in an undoped silicene nanoribbon. 
The band dispersion of silicene is sensitive to the applied 
electric field and posseses a topological phase transition 
from a topological insulator to a band insulator. Us¬ 
ing the lattice Green’s function technique in Lehmann’s 
representation we have derived different terms of spin 
susceptibility, which is related to the RKKY interaction. 
The RKKY interaction will have three different terms, 
namely, Heisenberg, Ising, and Dzyalosinkii-Moriya. As 
a result of the bulk-edge correspondence in the topolog¬ 
ical insulator phase, the gap closes at the boundary of 
the nanoribbon, and this shows its effect on the RKKY 
interaction for the impurities located on the edge. We 
have shown that, although the magnetic coupling decays 
fast and is short ranged, the RKKY interaction is sig¬ 
nificantly large in the case where two magnetic moments 
are located on the nanoribbon edges rather than sitting 
in the sheets. Remarkably, in the topological regime, the 
RKKY interaction is about 20 times greater than in the 
bulk insulator regime. 

One of the fascinating features in this case is the pos¬ 
sibility for the enhancement of the interaction by tuning 
the electric field, which opens an avenue to probe the 
interaction experimentally. In addition, we have demon¬ 
strated that the ground state of the system consisting of 
dilute magnetic moments on the edge can display a spiral. 
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ferromagnetic, or antiferromagnetic phase depending on 
the value of the out-of-plane electric field. This electri¬ 
cal tunability of the magnetic phases in dilute magnetic 
silicene can be investigated using current experimental 
techniques and can be of interest in the field of spintron- 
ics. 

As the final point, we want to address the experimen¬ 
tal accessibility of our proposals. The RKKY interaction 
is small such that measuring it experimentally is usually 
hard. However, using a high value of density of states on 
the edge of materials such as zigzag graphene nanorib¬ 
bons [13], or band tunability of materials such as bilayer 
graphene [16], increases the chance to measure this quan¬ 
tity. Zero gapped edge states in graphene only happen in 
perfect zigzag edges, however, there are both edge modes 
and band tunability in silicence. We should mention that, 
although we carried out the calculation for a zigzag sil¬ 
icene nanoribbon, since silicene is a Kane-Mele topolog¬ 
ical insulator [37], the edge modes are preserved even for 
rough edges. This fact increases the possibility that our 
theoretical predictions may be explored in experiment. 


Appendix A: Explicit form of the spin susceptibility 
of a nanoribbon 


In this section, we want to show the form of the spin 
susceptibility as a function of eigenvalues and eigenvec¬ 
tors of a ZSNR. Since we neglect the Rashba term, the 
Green’s function for spin-up and spin-down can be sepa- 

fg(') 0 \ 

rated in the matrix form ^ ( q (0 I ^t(l) 

are the spin components of G(r, r', s) and G(r', r, e), 
respectively. The various components of the spin suscep¬ 
tibility are given by r') = d6z/^’^(r, r'), 

where = Tr[cr^G^(r, r', 5)cr^G^(r', r, 5 )] and is 


given by 

= ^yy = gig't ^ gWi 

= g^g\ + g^g'i 

- g^g\) (Al) 

while the rest of the terms are zero. 

In order to calculate the spin susceptibility, we 
apply the Lehmann’s representation of the Green’s 

function (7). Each term of the multiplication 
of two Green functions can be rewritten in a 
compact form as = Sn,n'^nn'^nn' where 

Ki = C(r)V’r(r')C'(rOC''(r) and = 

[{s — s^-\- ir]){e — e^, -h Note that under the in¬ 

terchange of n and n', gsg',, —> En,n' Ki'^nn' and in 
a similar manner ^ 5 /^' converts to ^nn'* 

By applying these properties to the Eq. (Al) and using 
\im^^Q+{x^ if ])= P{l/x)^i7rS{x), and some straight 
forward calculations one can show 


X 




n^n' 


eh - e, 


+ X)€' (r)]. (A2) 

eti - e' 


In analogy to the the yields 


X^^(r,r') = 2 ^ M)_T(A)5Re[V.I(r)Cnr')V’I^(r')C'(r)] 


nj^n' 


t t 


+ 




On c. 


and finally in a similar manner, y^^ can be written as 




n^n' 


t t 

eh - et> 


(r')Cl' (r)]- 


t t 
eti - e\ 


(A4) 
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